Soft Elasticity in Nematic Liquid-Crystal Networks. 
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Liquid-crystal networks consist of weakly crosslinked polymers that are coupled to liquid-crystal 
molecules. The resultant hybrid system has rich elastic properties. We develop a phase field model to 
describe mechanical properties of a hexagonal liquid-crystal network. The hexagonal liquid-crystal 
network is found to have soft shear deformations. The elastic properties are predicted analytically 
and confirmed with numerical simulations. In addition the model naturally incorporates non-linear 
elasticity and dislocations or disclinations. 
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Soft elasticity is an exotic elastic phenomenon found 
in liquid-crystal networks which incorporates liquid- 
crystalline order and rubber elasticity [lj. The sim- 
plest liquid-crystal network (LCN) consists of weakly 
crosslinked polymers with the constituent rodlike liquid- 
crystal molecules in the nematic phase. A nematic 
LCN can relieve a shear stress via the additional de- 
grees of freedom that arise from the constituent nematic 
molecules. In addition to exhibiting rich elastic behavior, 
LCNs are being considered in new technological applica- 
tions that include artificial muscles Q, opto-mechanical 
actuators 0], and tunable mirror-less lasers 

Ideally, the nematic molecules are free to wiggle about 
their anisotropic axes. As a consequence, when a shear 
deformation is applied to a nematic LCN in a plane 
that involves the nematic director, the nematic molecules 
can tilt locally such that the system retains its equilib- 
rium configuration. This phenomena is referred to as 
soft elasticity Q and is in the same spirit as the pre- 
diction of Golubovic-Lubensky that an anisotropic glass 
that breaks spontaneously a continuous symmetry must 
have a vanishing shear modulus 5] . 

In this article we show that the phenomenon of soft 
elasticity is not limited to nematic gels [6] and smec- 
tic gels but can also be realized in networks with 
more translational order. Our system is distinct from the 
"soft crystal" phases of a liquid crystal fluid Q due to 
the topological constraints offered by the network. Our 
hybrid system thus shrinks or elongates spontaneously 
depending on the orientational distribution of the ne- 
matic molecules. Our model could be tested experimen- 
tally in a block copolymer system that can self-assemble 
in a liquid-crystalline environment into an ordered net- 
work. For example, a ABA triblock copolymer solution 
has been shown to physically crosslink when the ^4-type 
monomers are phobic to the liquid crystal solvent [§| . 

To capture the properties of a LCN we couple the net- 
work molecular shape to the orientational order of the 
background nematic molecules. The coupling is such 
that in the nematic phase, the network is elongated in 



a direction determined by the average orientation of the 
nematic molecules. In this state the network density field 
■0(r) is anisotropic. We assume that the network molec- 
ular lengthscale that is obtained from its characteristic 
wavenumber go is much larger than the dimension of the 
nematic molecules a. This assumption holds in the case 
of the ^45^4 triblock network system where the radius of 
gyration of the polymers R g = 1/qo >> a. This allows 
us to utilise the nematic director field 9(r) to describe 
the nematic molecules. 

Under these conditions we propose the phcnomcnolog- 
ical free energy functional 



F= J d 2 r{K{V8) 2 /2 + t^ 2 /2 



+ A^ 3 /3 + ^ 4 /4 + 



[q 2 iP + Ad xx iP + Bd vy ^ + Cd xy ^] 2 /2}, (1) 

where A{9) = cos 2 9 + n 2 sin 2 9, B{0) = sin 2 9 + k 2 cos 2 9, 
and C(9) = (k 2 — 1) sin 29 and where K is the Frank 
elastic constant for liquid crystals, k is the anisotropy 
ratio of the density fluctuations in the nematic state, r 
is a control parameter and A < is a phenomenological 
constant. In the limit k — 1, the network fluctuations are 
isotropic regardless of the nematic order and thus the two 
fields are uncoupled. 

The dynamics of the model are driven by the mini- 
mization of the free energy, i.e., 



dip(v,t)/dt = \J 2 [8F/5yb(r)}, 
d9(r,t)/dt = - [5F/S9(r)]+HTi(T,t), 



(2) 



where 77(1", t) is a noise field and fi is the intensity of 
the noise. For simplicity, the thermal fluctuations are 
only incorporated in the orientational field such that r\ is 
assumed to be a Gaussian random function, with a mean 
variance (r)(r, t), r)(r', t')) — 6(r—r')5(t — t'). The square 
of the intensity of the noise is linearly proportional to 
temperature. 

The equilibrium states for Eq. ([1]) can be determined 
in mean field theory by considering the minima of F. 
We can assume that the average molecular orientation is 



2 



along the y-axis [i.e., (9(r)) = 0]. In this instance the 
free energy functional simplifies to 

F m « J d 2 Y{{rld xx i: + r\d vv ^ + ql^) 2 /2 

+rV//2 + AV' 3 /3 + V' 4 /4}, (3) 

where r 2 = (B(6)), r\ = (A{9)) and the brackets imply 
an average over the liquid crystal orientational distribu- 
tion. Since the nematic director is assumed to be along 
the y— coordinate, from Eq. ([3]) we obtain that rii and r± 
define the periodicity of the network parallel and perpen- 
dicular to the nematic director respectively. 

The averaged free energy F m is in the form of the 
Landau-Brazovskii (LB) phenomenological theory [1(3] 
that describes phase transitions of a uniform system to 
a periodic state. The only difference is the scaling factor 
between the two coordinates. The mean field phase dia- 
gram of this model is well known ll|, [12[ . The approxi- 



mate equilibrium solutions of the density fluctuations are 
then similar to the crystalline solutions of the LB theory. 
For example the hexagonal phase is described by 

ip(x,y) = A t jcos VSq x/(2rx) cos [q y/(2r\\)] 

+ cos(g y/ni)/ 2 }, (4) 



where A t = 4(— A + \J — 15t + A 2 )/15 is the amplitude 
of the local density fluctuations. The sinusoidal solution 
Eq. ^ is the leading term approximation of a Fourier 
series and is only valid for > r >> — 1. 

At low temperatures, the nematic molecules have a 
preferred direction [6(r) =0, in our case] and thus we 
have r^ em = K and r" em = 1. In the high-temperature 
isotropic phase, the nematic molecules have no preferred 
direction, i.e., 0(r) is a random field. Thus we have 
r iso _ r iso _ y?( K 2 ZjZ l)/2. The change of r\\ and r±_ with 
temperature highlights elongations of the network lattice 
due to nematic ordering. An example of the network lat- 
tice distortion is presented in Fig[TJ Spontaneous elonga- 
tion of the network during a isotropic-nematic transition 
is a hallmark of LCNs Q. 

The continuum elastic free energy of our elongated 
hexagonal LCN is constructed from the lattice symmetry 
of the network and the coupling of its rotational defor- 
mation to the nematic director orientation. It is given 
by 

Eel = [Cn(u 2 xx + u 2 y ) + 2Ci 2 u xx u yy + AC 66 u 2 xy ] /2 



+ Di(9-uj xy ) /2 + D 2 {9 - Lo xy )u xy 
+ K{V9) 2 /2, 



(5) 



where — (dpUa + d a up)/2 is the symmetric strain 
tensor and co a p — {dpu a — d a u^)/2 is the antisymmetric 
strain tensor. u a /3 and uj a p describe respectively the rel- 
ative translations and rotations of the network due to an 
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FIG. 1: Equilibrium configuration of the local network density 
f/)(r) with anisotropy k = 2: (a) in the isotropic phase (fi = 
20) and (b) in the nematic phase (n = 1). The dark regions 
correspond to ip > and the light regions correspond to ip < 
0. (Inset: structure factor) 



applied deformation. The first line in Eq. ([5]) is the elas- 
tic free energy of a two dimensional hexagonal crystal 1 1311 
and the second line is the de Gennes elastic energy [l4j, 
which describes the coupling of the nematic molecules to 
local fluctuations of the network strands. 

We will determine the elastic moduli from our molec- 
ular free energy. We will assume that the distortions are 
small, so changes of the density amplitude ip(x,y) are 
negligible. We will choose the molecular axis and ne- 
matic director to be along the y-axis, where 9(r) = 0, 
and thus employ Eq. ([3]) as the effective free energy den- 
sity. A pure shear distortion is described by the den- 
sity ip — ip{x + Qy,y + C x )i where C is the elastic strain. 
In this state the shear free elastic energy per unit area 
(f shear = [F m (() - F m {( = 0)]/area) can be calculated 
by substituting a one mode approximation for ip [i.e., 
Eq. |gj)] into Eq. © to obtain 

f shear = ^ ( ^ + ^ q^C? + <D((*). (6) 

32 \r± r\\ J 

Eq. ([5]) can now be used to calculate Cqq recalling that 
for a pure shear u Qj a = 0, except for u xy — Q. This gives 
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Aht- (7) 



At low temperatures we have r^ em /V™ em = k, which 
means the shear modulus is strongly dependent on the 
network molecular anisotropy. Under a similar proce- 
dure, other physical elastic moduli such as the bulk 
modulus B u = 3v4(<7g/32 and the deviatoric modulus 
Cd = 3v4^<7 4 /16 are obtained. These do not depend 
on the molecular aspect ratio fl2T ] and note that when 
rii = rj_ then Cd — Cqq as expected for a hexagonal 
lattice. 

A pure rotational deformation is described by the 
strain tensor uj xy = £ and u a /3 — 0. The elastic en- 
ergy terms in Eq. ([5]) all vanish except the term with 
coefficient D\, In this state, tp = ip(x + Cy, V — The 
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evaluated free-energy density of the deformed state gives 

2 



D l = 



16 



rii 



A 2 n 4 



(8) 



From the difference between two simple shears along each 
coordinate we can determine Do.: 



Do 
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shear 



d( 2 
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U J she 
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The coupling moduli D\ and D2 vanish when the net- 
work density fluctuations are isotropic (where rj_ = rp), 
i.e, in the isotropic liquid-crystal phase or when k = 1. 
It is interesting to note that these calculations for the 
molecular moduli Cqq, D\ and D2 are consistent with 
the theoretical result of Warner and Terentjev |l[ for ne- 
matic elastomers derived from classical rubber elasticity. 

When the nematic director is allowed to relax to an op- 
timum state, then dE e i/dO = and thus Eq. (JSJ reduces 
to 

E fi = Cn (ul x + u 2 m ) /2 + Ci 2 u xx Uyy + 2C 66 ul y 
+ KD 2 2 (Vu xy f /(2Df) + K (Vcj xy f /2, (10) 



D 2 /(ADi) is the renormalized shear 



where Cqq 
modulus. 

From Eqs. ([7]), ©, and (0 we obtain the remarkable 
result that Cqq — 0. This means the nematic direc- 
tor relaxes to cancel out the elastic energy cost for a 
shear deformation. In this limit, the nematic molecules 
"wiggle" until the deformed structure is compatible with 
the boundary conditions. The stability of the crystalline 
state at zero shear modulus requires higher orders terms, 
like 0(£ 4 ), in the elastic free energy. 

To confirm these approximate analytic calculations, we 
numerically solve Eq J2]). To examine soft elasticity, we 
impose a steady shear deformation by adding an advec- 
tive term on the dynamics, i.e., d/dt — > d/dt + V x d x , 
where the velocity, V x = 7?/, has a gradient in the y-axis. 
We set the initial nematic director orientation such that 
it is along the shear flow direction and we also shift our 
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FIG. 3: The free-energy density evolution after commence- 
ment of a shear flow with a shear rate 7 = 0.002. The system 
is shown at different molecular ratio: k = 1.0, 1.35, and 1.7, 
for the dotted, dashed, and solid lines respectively. (fi = 1.0) 



periodic boundaries to be consistent to the shear flow 
using the Lees-Edwards method [3]. We choose the fol- 
lowing numerical parameters: r = —0.03, A = —0.9, and 
K = 1, such that the network forms a hexagonal lattice. 

The influence of the shear on the total free energy den- 
sity F s is displayed in Fig. [2] for three values of the ther- 
mal noise strength, fi. At the lowest temperature the free 
energy is flat at small strains up to 'yt sa 0.25. Thus the 
system has the same free energy as the undeformed state 
at small strains. However at higher temperatures, F s 
increases for all shears as it is expected for a conventional 
network. These results are consistent with Eq. ([TO)) , since 
the coupling constants D\ and D2 vanish with increasing 
temperature. Thus soft elasticity vanishes with increas- 
ing temperature. We also examined the influence of the 
network anisotropy on the shear deformation and the re- 
sults are displayed in Fig. [3] When k = 1 the two fields 
are decoupled, and the elastic response is that of classi- 
cal solid. When k = 1.7 a soft regime is observed up to 
yt s» 0.3 as indicated by the flat free-energy density. 

Finally it is interesting to consider the configuration of 
the network density V'( r ) an( i the spatial average orien- 
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a static strain. For small strains (0(r)) relaxes to the 
value determined by Eq. (fTTj) and the total free energy 
decays to the value of the undeformed state in agreement 
with our mean field predictions. We tested our numerical 
computations for —0.03 < r < —0.3 where the hexagonal 
phase is stable, and these simulations indicate that the 
analytic results are exact as r — ► 0_, with small correc- 
tions at larger |r|. The limit r — > 0_ is also where our 
sinusoidal density approximation is valid. 

In summary, a model of liquid crystal networks was 
presented and shown by analytical and numerical meth- 
ods to reproduce soft elasticity as a function of temper- 
ature and molecular shape. In addition the numerical 
simulations provide evidence of unusual non-linear yield- 
ing mechanisms which provide avenues of future research. 
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FIG. 4: Evolution of the free-energy density F s and the ne- 
matic director orientation after commencement of a steady 
flow with shear a strain rate 7 = 0.0025. The snapshots of 
the density configuration correspond to the data at the given 
strains. The dashed line corresponds to the linear elastic the- 
ory. (k = 2,h= 1.0) 



tation of the nematic molecules (0(r)) as the system is 
being sheared. The snapshots in Fig.[|]show that the net- 
work structure changes with increasing shear strain even 
in the elastic soft regime. A uniaxial hexagonal lattice 
has a continuous set of structures which within our 
mean field theory have the same free energy. The switch- 
ing of (f5(r)) is different from the rotational component 
of the shear which is given by oj xy = jt/2. The initial 
dependence of (0(r)) on the shear can be calculated by 
minimizing the elastic energy E e i for a simple shear de- 
formation applied along the nematic director. This gives 

6 soft (t) = (£>! + D 2 )C(t)/(2Di) = C(*)/(l - « 2 ) (11) 

in the long wavelength limit. As shown in Fig. 2] this pre- 
diction works quite well. After this soft regime the strain 
energy increases until a yield occurs, e.g., for k — 2 at 
jt ~ 0.75. Interestingly this first yield occurs without 
the nucleation of mobile dislocations as would occur in 
a normal crystalline material. Of course at higher strain 
dislocations do eventually appear as can be seen in Fig. |4] 
We have also made numerical calculations in the case of 
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